Absence of static orbital current magnetism at the apical oxygen site in HgBa 2 Cu04 +< 5 

from NMR 
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The simple structure of HgBa2Cu04+i (Hgl201) is ideal among cuprates for study of the pseudo- 
gap phase as a broken symmetry state. We have performed 17 nuclear magnetic resonance (NMR) 
on an underdoped Hgl201 crystal with transition temperature of 74 K to look for circulating orbital 
currents proposed theoretically and inferred from neutron scattering. The narrow spectra preclude 
static local fields in the pseudogap phase at the apical site, suggesting that the moments observed 
with neutrons are fluctuating. The NMR frequency shifts are consistent with a dipolar field from 
the Cu +2 site. 

PACS numbers: 74.25.nj, 74.72.Kf, 76.60.Jx 
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A distinguishing feature of cupratc superconductors, 
other than their high transition temperature, is the ex- 
istence of a normal state pseudogap in underdoped ma- 
terials. However, the origin of the pseudogap is not yet 
understood. There are many suggestions that it is con- 
sistent with a broken symmetry state such as staggered 
fluctuating currents[l , orbital currents which conserve 
translational symmetry 2 4 , <i-density wave [5 , and var- 
ious other broken symmetries [7]. 

Polarized neutron scattering in YBa9Cu^O7_A[8lU0|. 
HgBa 2 Cu0 4+I 5 (Hgl201)[ffllI2], Bi 2 Sr 2 CaCu 2 8+(5 [I3], 

and La 2 _ 2 Sr :E Cu04 (short ranged) [14] show the appear- 
ance of broken time reversal symmetry that correlates 
with the onset of the pseudogap phase with a magnetic 
moment tilted away from the crystalline c-axis. A muon 
spin relaxation (i/SR) experiment in zero field was un- 
able to detect this magnetism [15] although screening ef- 
fects may account for muon insensitivity.[16 To account 
for neutron data, extensions of 2D orbital currents to 
3D models have been proposed which include the apical 
oxygen[17 as well as a quantum mechanical calculation 
[18j which has an out-of-plane component with a tilted 
local field. In this context, magnetic fields at the apical 
oxygen site, Fig. 1(a), up to 200 G have been calculated 
which should be measurable by NMR. [TO] 

Both NMR and nuclear quadrupolar resonance (NQR) 
experiments have previously attempted measurement of 
orbital currents: 89 Y NMR in Y 2 Ba 4 Cu 7 Oi5 +( 5[20] and 
135 Ba NQR in YBa 2 Cu 4 8 .[2T] However, 89 Y is in a sym- 
metry position which is insensitive to orbital current or- 
dering and neutron data is not available in these mate- 
rials for direct comparison. We report here our investi- 
gation of local fields in Hgl201 with 17 NMR for which 
the apical oxygen is sensitive to orbital currents and po- 
larized neutron scattering show evidence for magnetic or- 
dering. 



Until recently, NMR results for Hgl201 have been 
constrained to c-axis magnetically aligned powdered 
samples, [22 -25] which are more easily annealed and have 
a large NMR signal. However, these powder samples are 
disordered in the afe-plane and they can have a larger 
spectral linewidth due to imperfect alignment. In order 
to investigate possible effects of orbital currents, high 
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FIG. 1: (a) Orbital current ordering (green arrows) in the 
CuOe octahedron in Hgl201 proposed [171119] as explanation 
for neutron scattering results, (b) The profile of the currents 
and the effective magnetic fields (blue symbols) produced. 
The apical oxygens will experience different field depending 
on lattice position with 0(2)a above and 0(2)b below the 
CUO2 plane while the local fields will cancel at the O(l) site 
by symmetry, (c) & (d) A view of the CuO-plane along the 
c-axis, showing the local fields (blue arrows) at the apical 
oxygen (above and below the Cu nucleus) at two different 
orientations: Hq II a and Ho — 45° between the a and 6-axes. 



quality single crystals are necessary. These samples have 
become available[26j [27] and the results of 199 Hg and 
63 Cu NMR[28] [29] have been reported. We isotope ex- 
changed 17 and performed high resolution NMR fre- 
quency shift measurements for both the planar 0(1) and 
apical 0(2) oxygen, with narrow linewidths of approxi- 
mately 12 kHz and 5kHz respectively at Ho = 6.4 T, for 
a moderately underdoped sample with superconducting 
transition temperature of 74 K, as compared to 95 K at 
optimal doping. 

If the pseudogap state breaks time reversal symmetry 
from orbital currents, as shown in Fig. 1(a), the two api- 
cal oxygens in each unit cell will experience oppositely 
directed field components, Fig. 1(b). Furthermore, the 
orientation of these magnetic fields are locked at 45° from 
the a or 6-axes. These field components have an angular 
dependence with a maximum value, H max , between the 
a and 5-axes, H max /\/2 along each axis, and zero parallel 
to the c-axis. With applied magnetic field in the a6-plane 
the orbital currents will produce a splitting of the cen- 
tral transition with a period of 7r for rotation about the 
c-axis. 

Single crystals were grown by self flux method [221 [STJ 
then heat treated for 17 exchange and annealing. A 20 
mg sample was held at 600° C for 12 hours in 1 atm of 
90% 17 and then, in a continuous step, lowered to 300°C 
for 10 days. The annealing resulted in an underdoped 
superconducting transition temperature, T c , of 74.0 K 
with a transition width, AT C , of 3 K, compared to 74.5 
K and 3 K, respectively, before exchange. 

Nuclear magnetic resonance was performed in fields of 
Hq = 6.4 and 16.5 T and temperatures from 4.2 to 295 
K. Both stoichiometric sites were readily observed and 
the resonant frequency was taken from the peak of the 
spectra. Spectra were measured by a standard Hahn echo 
sequence, 7r/2 — r — it, with a typical ir/2 pulse length 
of 2 /is, sufficient to excite both sites. For full spectra, 
including satellite transitions, Fig. 2, a frequency sweep 
method was used. From the full spectra the quadrupolar 
splitting was found to be vq = 1.08 MHz for both sites 
with the axis of symmetry along the a and c-axes for 0(1) 
and 0(2) respectively. The 0(2) site exhibits axial sym- 
metry while the 0(1) site has less than axial symmetry 
with an asymmetry parameter rj = 0.40. With Hq _L c, 
we were able to separate the 0(1) and 0(2) resonances 
taking advantage of their rather different spin-spin relax- 
ation times, i.e. 0(2) has a much longer relaxation time 
and consequently, for sufficiently long r, 0(1) loses phase 
coherence, leaving only an 0(2) contribution to the signal 
amplitude. 

The NMR spectrum frequency shift v — jHo = 
K(T, Ho) r yHo, where v is the resonant frequency, 7 is 
the gyromagnetic ratio, and K(T,H ) can be written as 
a sum of components such that, 

K(T,H ) = K s (T)+K o +K q +Att(1-D)M(T)/H , (1) 
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FIG. 2: NMR spectra at H = 16.5 T and T = 100 K, includ- 
ing satellite transitions, for Hq || c and Ho \\a. The transitions 
are labeled na where n — 1,2 for 0(1), 0(2) and a — A, B, 
C for the central, first quadrupolar, and second quadrupolar 
transitions. The central transition and four satellites can be 
seen for three distinct oxygen sites, apical, planar with Cu-O 
bond along the field, and planar with the Cu-O bond per- 
pendicular to the field (primed). The planar sites should be 
degenerate with Ho 1 1 c however the crystal is misaligned by 
«10 degrees, evident by two sets of satellites in this orienta- 
tion. The higher order satellites 2C and 1C are omitted for 
clarity. 



and K S (T) is the spin shift, called the Knight shift; K 
is the orbital shift; Kq is the quadrupolar shift; D is the 
demagnetization factor (dependent on sample geometry) , 
and M(T) is the magnetization due to superconduct- 
ing screening currents, all of which can be anisotropic, 
although only K S {T) and M(T) have temperature de- 
pendence. As is usually considered for cuprates. pffl I3"2"] 
we adopt a single-component spin susceptibility model, 
K s = Axs(T), where A is the transferred hyperfinc field 
for a particular nucleus and Xs(T) is the electronic spin 
susceptibility. 

Our main results are the orientation and tempera- 
ture dependence of the 0(1) and 0(2) frequency shifts 
of the central transitions, Fig. 3. For all orientations, 
the 0(1) shift exhibits the classic manifestation of the 
pseudogap [3"5I 13"4"] . decreasing from above room temper- 
ature down to our minimum temperature T = 4K with 
no clear indication of the transition to the superconduct- 
ing state. In our sample, onset of the pseudogap phase 
can be estimated from polarized neutron scattcring^ll 
to be at T w 270 K. The shift of 0(2) behaves signifi- 
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FIG. 3: Temperature dependence of the planar 0(1) and api- 
cal 0(2) oxygen shifts with Ho || c, Ho _l_ c with Ho || a and 
with Ho at 45° to the a and 6-axes. For Ho || a, there are two 
inequivalent O(l) sites. [30] The O(l) shift shows the classic 
decrease of K s characteristic of the pseudogap. The 0(2) shift 
has a minimal K s component which allows for clear observa- 
tion of the reversible diamagnetism in the superconducting 
state. 



cantly differently. Above T c there is little temperature 
dependence, while below T c there is a sharp onset for 
a decreasing temperature dependent shift which we at- 
tribute to superconducting diamagnetism. 

We look for the existence of orbital current order, con- 
tributing to K , in two different ways. First, we compare 
the field dependence of the shift along the c-axis with that 
in the afr-plane along the a-axis. Secondly, we measure 
the angular dependence of the shift and the linewidth for 
a magnetic field within the a6-plane. 

The spin susceptibility is a strong function of tempera- 
ture, evident from the 0(1) shift, decreasing by over half 
its value from room temperature down to T Cl Fig|3] Thus, 
the small temperature dependence of the 0(2) shift in the 
normal state indicates that the transferred hyperfine field 
from the CuO-plane is correspondingly small. Therefore, 
the major contributions to the shift of the apical oxygen 
can be taken to be, 



(2) 



K (2)(0,H ) = K o (0)+K y Q z> (9,H o ) 



(2) 



where 9 is the angle between the magnetic field and 
the c-axis. The central transition is unaffected by the 
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FIG. 4: Central transitions at T = 100 K at 15° intervals from 
6 = to 90°, or Ho \\ c to Ho || a respectively at Ho — 16.5 
T. For the planar oxygen, the angular dependence shows 
the rhombic symmetry of the orbital, spin, and second or- 
der quadrupolar shift. The apical oxygen has axial symmetry 
along the c-axis. Inset: Angular dependence of the 0(2) sec- 
ond order quadrupolar and orbital shifts, Eq.2, given by the 
blue curve. 



(2) 

quadrupolar interaction to first order, so Kq is the 
second order shift given by the blue curve in the inset 
to Fig-El v = v 2 Q /2v (l - 9cos 2 0)(3cos 2 e - 1), where 
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FIG. 5: Frequency shifts of the apical oxygen in the under- 
doped sample with orientations Ho 1 1 c and Ho 1 1 a as a func- 
tion of temperature relative to the value at T c . The shift be- 
low T c can be attributed to superconducting diamagnetism. 
Inset: Above T c the shifts are deconvolved into isotropic and 
axial components showing that only AK ax is temperature de- 
pendent. 



vq is the Lamor frequency and vq is the quadrupolar 
splitting from the satellite spectra, Fig(2] Our data 
for 0(2) closely follow this behavior, which after sub- 
traction leaves the orbital components. Circulating or- 
bital currents would shift the resonance frequency an 
amount independent of applied field. At T = 100 K, 
and H a = 6.4 (16.5) T, we find the difference AK = 
K o (0 = 90°) - K {6 = 0°) = 0.025% (0.027)%, to be 
almost field independent placing a limit on a possible ex- 
cess field from orbital currents of 2.1 G in the a6-plane. 

For an independent estimate of the effect of orbital 
currents in the a6-plane, we have measured the angular 
dependence of the shifts at a magnetic field of Hq — 6.4 T 
rotated about the c-axis. As can be seen from Fig. 3(b), 
the apical oxygen NMR frequency has no in-plane an- 
gular dependence within our resolution of 2 G, indi- 
cating that there is little, or no, additional contribu- 
tion from orbital currents. A more stringent test, if 
the orbital currents have the symmetry shown in Figjl] 
is to examine the NMR linewidth for broadening that 
might reflect a small frequency splitting. Our measure- 
ments at all temperatures of the angular variation of the 
linewidth is less than 0.5 kHz which places a bound on 
orbital current fields of 0.4 G. According to Lederer and 
Kivelson.|19j orbital currents inferred from neutron scat- 
tering results [HI HQl |TT] should lead to local fields of order 
H m ax ~ 200 G oriented as indicated in Figjl] Taken to- 
gether, our measurements of the angular variations of 
the apical orbital shift and NMR linewidth, in-plane and 
out of-plane, place substantial constraints on predictions 
from the theory for a static orbital field to be less than 
0.4 G. 

Examining the 0(2) shift in greater detail reveals that 
the 0(2) nucleus has a small temperature dependent spin 
shift in the normal state such that from room tempera- 
ture to T c , AK S equals (-0.0035,-0.0035,0.0065)% in 
the crystal basis (a,b,c), FigJ5J The temperature de- 
pendence of the shifts can be deconvolved into isotropic, 
Ki SO , and axial, K ax , components where the latter has 
the angular dependence of the second order Legendre 
polynomial. For 0(2), the shift is axially symmetric 
around the c-axis and the shifts are given by, 



(K c + 2K ab )/3 



K, 



K ax = (K c - K ab )/3 



(3) 



This decomposition shows that all of the temperature 
dependence resides in the axial component of the spin 
shift, Fig [5] The transferred hyperfine field from the 
Cu +2 ion and contributions to the ligand 0~ hyperfine 
field, can come from several sources 35 : 



A c + A. 
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A c is the Fermi contact term, A cp is from core polariza- 
tion, A p is from unpaired p electrons (only p is relevant 
for the 0(2) nucleus, although d and / would also be 
axial), and A dip is the direct dipolar coupling from spin 
on the Cu +2 ion. From the temperature independence of 
Ki SO we infer that the hyperfine fields A c and A cp are 
negligibly small, corresponding to temperature depen- 
dent shifts less than 3xl0~ 4 %. Core polarization comes 
from the unpaired p electron and A p ex A cp . The ratio 
can be estimated as A cp /A p = (10/4)re where re is an 
empirical constant. For atomic fluorine (the same elec- 
tronic configuration as O - ), re = —0.1 36J, giving a ratio 
of A cp /Ap — 0.25. This would place an upper bound 
on A p corresponding to at most a 40% contribution to 
AK ax , suggesting that the majority of its contribution 
comes from a dipolar field such as a magnetic moment 
on the Cu site. [HJ |37j 

Presuming the dipolar effect at the 0(2) site also ex- 
ists at the planar 0(1) site, and taking the 0(2) shift 
from 300 K to T c , AK aXt0{2) = AK dip>0(2) = 0.0035%, 
we find for the 0(1) site AK dip0 ^ = 0.021%, noting 
that it has two Cu nearest neighbors and using known 
atom distances. [35] This shift agrees quite well with our 
measured anisotropic shift of 0(1), AK ax ^ (i) — 0.018%, 
indicating that K ax Q(i) has a significant dipolar contri- 
bution from the Cu +2 site as well. 

The hyperfine fields at the apical oxygen site are suf- 
ficiently small that we can readily observe the shift in 
the superconducting state from superconducting diamag- 
netism, Fig(5] Taking into consideration the demagneti- 
zation factor that is dependent on the overall crystal di- 
mensions, we determine the superconducting penetration 
depths, A Q , where a = || or _L denotes the orientation of 
the field relative to the c-axis, AK = K(T C ) - K(T = 
0) = 47r(l — D a )M a /Ho. The crystal dimensions give 
demagnetization factors D± = 0.67 and Dm = 0.165 |39j . 
The magnetization, M a (H), can be expressed by: 



M a (T) 



$0 l n Hc2 

87rA?,(T) H ' 

$0 1„ H_c2 

87rA,|(T)Ax(T) m H 
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(5) 



where Ai SO (A ax ) are the isotropic (axial) hyperfine fields, 



where $0 is the magnetic flux quantum and H C 2 is the up- 
per critical field taken to be 80 T. [4DHI2"] Solving for the 
penetration depths we find Am (T = 0) = 1, 760±60 A and 
X±(T = 0) = 1.7 ± 0.2/xm. Earlier results in near-zero 
applied field from imaged Josephson vortices ^43 found 
Aj_(T = 0) ~ 8/zm , significantly larger than ours. How- 
ever, A|| and Aj_ from susceptibility on p owders [4i" 1 142 | [44] 
are consistent with our work. This comparison shows 
that there is no significant variation in the penetration 
depth for fields up to Hq = 6.4 T indicating the minimal 
role of non-local effects due to Doppler contributions, also 
referred to as the Nonlinear Meissner Effect. [45] 

Absence of magnetic field shifts from 0(2) NMR rules 
out the existence of static circulating orbital currents 
at the apical oxygen site in HgBa2Cu04-|-,5. However, 



discrepancy between detection of magnetic ordering by 
neutrons which are absent in NMR measurements is 
not unique to Hgl201. For example antiferromagnetic 
ordering has been observed by neutron scattering in 
UPt3 [IB] but not detected in NMR measurements which 
was suggested to be a consequence of the effects of 
thermal fluctuations and different time scales for these 
experiments. [47| 148] Roughly speaking these are 10~ n s 
for neutrons, [5] 10" 8 s for ^SR,[T5] and 1CT 6 for NMR fre- 
quency shifts. In the context of Hgl201 this was noted 
by Li et al. [H] 

In summary, we have measured the temperature, field, 
and angular dependence of the narrow planar and apical 
17 NMR spectra in an underdoped Hgl201 single crys- 
tal. There are no indications of static magnetic fields 
from circulating orbital currents proposed theoretically, 
within our accuracy of 0.4 G. The apical oxygen spectra 
have a small temperature dependence above T c which 
has axial symmetry, consistent with a direct dipolar field 
from the Cu +2 site. 
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